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We have analyzed the transmission properties of pulses through one-dimensional periodic 
structures in order to systematically explore the best conditions to achieve the maximum 
delay with the minimum possible distortion. In the absence of absorption and no layer 
variation, the transmission coefficient t m can be well approximated by a sum of Lorentzian 
resonances. The ratio between their width Tr and their separation Acjr is a crucial parameter 
to characterize the distortion of the transmitted pulse. For typical values of the parameters 
used in telecommunications and high index of refraction contrasts n2/ni, the distortion of 
the transmitted pulse is unacceptably large for frequencies near the edge of the transmission 
window. We estimate fractional delays achievables in terms of the central frequency used 
and the pulse bandwidth. 


I. INTRODUCTION 

Slow light phenomena has reeently attraeted a great deal of interest. Many researehers have 
delayed optieal pulses in a great variety of physieal systems sueh as atomie gases [[I]-l3l, optieal 
fibers lU-O or photonie-erystal deviees IITl-fTOll. 

The ability of delaying optieal pulses has many potential applieations in the field of high-speed 
optieal proeessing sueh as random-aeeess memory, data synehronization or pattern eorrelation. 
In these communications applications, where large bandwidth systems are required in order to 
respond fast to the short pulses which transport the information, the pulse of light must be delayed 
several times the pulse duration without significant distortion. 

Recently, Mok et al. IfTTI managed to delay optical pulses through the use of solitons launched 
near the band gap edge of a fiber Bragg grating (FBG). They demonstrated the propagation of 
0.68 ns optical pulses travelling at Vg = 0.23/c in a 10 cm FBG leading to delays corresponding 
to about 2.4 pulse widths. Sarkar et al. ifT^ obtained a fractional pulse delay of 200 % for 8 
ps pulses in a GaAs/AlGaAs multiple quantum well sample and Novikova et al. [TTSl managed 
to slow light in Rb vapour with minimal loss and fractional pulse delays greater than 10. In a 
numerical simulation of light propagation, Boyd et al. IIT4l obtained a delay of 75 pulse lengths 
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under realistic laboratory conditions. 

Our group has studied the traversal time of electronic and photonic pulses concentrating on 
the case of exponentially decaying wave functions, i.e., tunneling experiments in electronic and 
frequency gaps in optical transmission [fTSll . Our previous aim was to study how fast a wave could 
travel in a given region. For that purpose, we developed a method to calculate the traversal time 
including finite size effects, i.e., taking into account the specific form of the pulse. We found that 
the traversal time r for an incident wave packet with Fourier components $ 1 ( 0 ;) which traverses a 
structure of length L can be written as IfTSlI 


du \t{u)nUu)\^ T,{U) 

du 


( 1 ) 


where \t{uj)\‘^ is the transmission coefficient and ti{uj) is the so called Biittiker-Landauer time or 
phase time, which corresponds to the frequency derivative of the phase of the complex transmis¬ 
sion amplitude ifT^ 


n{uj) = 


diptjuj) 

doj 


( 2 ) 


In this paper we analyze under which conditions we can slow down light pulses as much 
as possible (large fractional delays) with minimum distortion in one-dimensional (ID) periodic 
structures, such as FBGs, usually employed in optical communication systems. The most impor¬ 
tant parameter to quantify how slow light is traversing a region is the fractional delay FD of our 
transmitted pulse, which is defined as 


FD = 



( 3 ) 


where At is the width of the incident pulse. We have used Eq. Q to calculate the fractional delay 
of the transmitted pulse, and we have found that under most circumstances the classical calculation 
constitutes an excellent approximation. Within this approach the delay time is given by 


FDel = —- 

/At \V„ C 


( 4 ) 


where Vg is the group velocity of the transmitted pulse. 

As the group velocity tends to zero as we approach the edge of the transmission window pro¬ 
duced by a periodic structure, it is natural to use frequencies near this edge to slow down light as 
much as possible. However, we have found that, for high index of refraction contrasts n 2 /ni, the 
pulse distortion is unacceptably large near the edge and usually does not compensate to work in 
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this region. For a given degree of distortion and high eontrasts, it is better to use short pulses near 
the eenter of the frequeney window to obtain large values of the FD. 

We employ the eharaeteristie determinant method [fTTII to ealeulate the transmission eoeffieient 
through a layered periodie strueture [fTSl . This method eonstitutes an exeellent tool to derive 
analytieal expressions and to obtain numerieal results for all the relevant quantities in periodie 
arrangements. We have eonsidered an idealized case with no absorption and no layer variation. 
The inclusion of any of these effects is difficult and outside the scope of this work, although it will 
be the subject of future work. In the discussions, we set the limits of validity of our present model 
and discuss possible effects due to absorption and layer variability. 

The plan of the work is as follows. In Sec. II we present the ID periodic structure used in 
our numerical calculations and describe in detail its transmission coefficient obtained with the 
characteristic determinant method. In Sec. Ill we study the pulse propagation through the ID 
grating depending on the position of its central frequency in the transmission window, and in Sec. 
IV we calculate numerically the transmission coefficient and the fractional delay of our pulse as 
a function of the frequency bandwidth. We also analyze the convenience of using central pulse 
frequencies near the center or the edge of the transmission window. Finally, we summarize our 
results in Sec. V and present an Appendix with details of the calculation of the transmitted pulse 
through a periodic arrangement of Lorentzian resonances. 


II. TRANSMISSION COEFFICIENT ANALYSIS OF THE ID PERIODIC STRUCTURE 

The slow light structure considered throughout the paper is a periodic arrangement of layers 
with index of refraction ni and thickness Li alternating with layers of index of refraction n 2 and 
thickness L 2 . The wavenumbers in layers of both types are fcj = ojriijc for z = 1,2 and the 
grating period is a = Li + L 2 . We have chosen niLi = 712 ^ 2 . The periodic modulation of the 
index of refraction can be constructed in the core of a short segment of optical fiber constituting 
a fiber Bragg grating. Such structures present a photonic band gap, where light pulses are totally 
reflected, and are typically used as inline optical filters (see Fig.[^. In our numerical calculations 
the FBG consists of 40000 Si layers of length 158.8 nm and index of refraction 3.48 alternating 
with 39999 air layers of 553.2 nm wide. The grating period is then 712.0 nm and the total size of 
our arrangement 28.5 mm. This corresponds to a realistic structure used in slow light experiments 
Omi, but the results are of general validity. 
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Figure 1: A fiber Bragg grating with an uniform periodic modulation of the index of refraction. 


To obtain the transmission coefficient of the previous structure |tArp, we employ the character¬ 
istic determinant method IfTTI . This formalism allows one to express the transmission coefficient 
of a wave propagating in a ID periodic structure through the determinant D^, which depends on 
the amplitudes of transmission and reflection of a single cell. Via this determinant method, the 
inverse of the transmission coefficient for this structure can be written as IfT^ 


I^Afl ^ — 1 + 


^2 _ k2 

2kik2 


sin(/c2-f^2) 


sm\NPa/2) 

sin^(/Sa) 


( 5 ) 


where iV — 1 is the total number of layers and (3 plays the role of quasimomentum of the system, 
and is defined by 


//^2 _|_ ^ 2 \ 

cos(/?a) = cos(fciLi) cos(/c 2T2) - ^ sin(fciLi) sin(/c2-f^2) 

V 2fciK2 / 


( 6 ) 


When the modulus of the RHS of Eq. ([^ is greater than 1, (3 has to be taken as imaginary. This 
situation corresponds to a forbidden frequency window. 

In Fig.j^a) we have represented the transmission coefficient of our grating as a function 
of the frequency oo for N = 10 to illustrate its behavior. The realistic grating corresponds to Fig. 
I^b) where N = 80000 (with 28.5 mm length). There are so many peaks that they appear as a 
black region in the figure, and they cannot be resolved in the scale shown. In order to study 
in detail as a function of the position in the frequency transmission window we plot on a much 
larger scale a region near the center. Fig. [^a), and a region near the edge. Fig. [^b). The center of 
these two regions, cuoi = 1754 THz and a;o 2 = 2225 THz, are marked with arrows in Fig.|^b) and 
they correspond to wavelengths within the infrared spectrum (Aqi = 1073.6 nm and A 02 = 846.3 
nm) for the second transmission window. In percentage terms, the frequencies cuoi and 0202 are at 
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Figure 2: Transmission coefficient |fArp as a function of the frequency a; for (a) = 10 and (b) N = 

80000. The regions around the frequencies tuoi = 1754 THz and 0102 = 2225 THz marked with arrows are 
expanded in Fig. 

45 % and 1 % of the total window width from the upper frequeney window edge, respeetively. 
Similar results will be obtained for the other frequeney windows. 

The dashed lines in Fig. [^a) and |^b) eorrespond to Lorentzian fits of the eentral resonance 
for each case. We note that the Lorentzian curves fit fairly well the peaks in the transmission 
coefficients in all cases studied. We can in fact justify this by having a close look into Eq. Q 
where the term — 1 is the product of a highly oscillating function, sm‘^{Nf3a/2), and factors 
independent of the number of layers and which vary very slowly on the scale of one peak. The 
oscillating function can be approximated around each of its minimum by a quadratic term, result¬ 
ing in a Lorentzian curve for the transmission coefficient. We will use this expansion later on to 
calculate the width of each peak. 

One observes that the resonances are narrower and closer to each other near the transmission 
window edge. From now on, we will denote by Tr the energy width of each individual Lorentzian 
resonance and by AcUr the frequency separation between two consecutive resonances. The ratio 




























6 


CM 







0.8 


0.6 


0.4 


0.2 


. 0.0 


1.0 


1754.512 1754.528 1754.544 


2225.875 2225.880 2225.885 


Q) (THz) 


Figure 3: Transmission coefficient versus w for a region near the center (a) and the edge (b) of the 
transmission window plotted in Fig.j^b). The dashed curves are Lorentzian fits for a single peak. 

Fr/^AcOr, which we will see that eonstitutes an important parameter of the problem, deereases as 
we approaeh the window edge. For this reason, the overlap between Lorentzians and the transmis¬ 
sion in the midpoint between peaks is very small in this region. 

We ean now obtain analytieal expressions for the width Fr and the separation AcOr. The term 
within braekets in Eq. ([^ only depends on the properties of one layer, while the quotient of the sine 
funetions eontains information about the interferenee between different layers. The transmission 
eoefficient is equal to 1 when sin {Nf3a/2) = 0 and (5 is different from 0. This condition 


occurs for 



77 = 1,...,A/2-1, 


(7) 


and it corresponds to a resonant frequeney. Considering that eaeh resonanee is a Lorentzian fune- 
tion eentered about the resonant frequency ccr with full width V^/h, we ean write the inverse of the 
transmission eoeffieient for one resonanee as 



( 8 ) 


Expanding the sine in the numerator of Eq. Q near a resonanee we have 



(9) 


so, Eq. Q ean be rewritten as 



( 10 ) 
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where we have taken into aeeount the definition of the group veloeity Vg in terms of the quasimo¬ 
mentum (3 


= 


6u 

Ip' 


( 11 ) 


We ean derive an expression for the frequeney distanee between two eonseeutive resonanees 
AcUr by taking finite differenees in Eq. (|7]) (resonanee eondition) 


27r 

AcUr = V„A3 = Ug-—, 
® ^Na 


( 12 ) 


and after eomparing Eqs. @ and ([T^ we ean identify 


'2\ sin(/)a)| 


hAUr 


TT 


_ u2 

2kik2 


sin(/c2E2) 


-1 


(13) 


As we will see, the ratio Tj-fhAur is a very signifieant parameter whieh is independent of the 
number of layers of our grating. In Eig. we plot this ratio, ealeulated via Eq. ( [T3| ), as a funetion 
of the frequeney oo (left axis, squares for A = 10 and solid line for N = 80000). One observes 
that in the eenter of the transmission windows the ratio Tj./hAu^ reaehes its maximum, whieh is 
roughly 1.0 in our ease. This ratio also tends to zero as we approaeh the optieal mirror region, 
what has important eonsequenees on the distortion of the transmitted pulses, as we will see. 

It is well known that the group veloeity tends to zero when approaehing the edge of the trans¬ 
mission windows. This is the reason why most authors have tried to work in this region to slow 
down light as mueh as possible. In Eig.|^we also show the ratio of the group veloeity to the speed 
of light in vaeuum Vg/c (right axis, dashed line) versus the frequeney u. As expeeted, this ratio 
tends to zero as we approaeh the transmission window edge. 


III. PULSE PROPAGATION THROUGH THE TRANSMISSION WINDOW 

Eet us now study how an ineident gaussian pulse \I/inc(^) traverses our ID grating depending 
on the position of its eentral frequeney in the transmission window. We will eonsider the two 
frequeneiesdesoribedinEig.[^b),a;oi = 1754 THz (near the eenter) and a;o 2 = 2225 THz (near the 
edge). In both oases, we have ohosen a pulse width of 68.4 ps whieh oorresponds to a bandwidth 
of about 41.3 GHz. 

Eor these broadband pulses we ean approximate the transmission amplitude of our grating by a 
periodio funetion of u that ean be written as the sum of a given funetion displaeed periodioally by 
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Figure 4: The ratio Fr/^Awr (left axis, squares for N = 10 and solid line for N = 80000) and the fractional 
group velocity Vg/c (right axis, dashed line) versus the frequency uj. 


a fixed amount 

(14) 

n 

As before, Au^ is the frequency separation between consecutive resonances. As we saw in the 
previous section, we can consider, to a very good approximation, / to be a Lorentzian function 
of width Fr. Under these conditions, we have shown in the Appendix that the transmitted pulse 
consists of an exponentially decreasing series of similar pulses given by 




2n 




t 


27rn\ 

AUr) 


(15) 


where (7 is a normalization constant. If we neglect overlap between consecutive peaks, it is easy 
to deduce from the previous expression that the ratio r between the intensity of two of these peaks 
is equal to 


r = exp 


27rrr 


hAUr 


(16) 


According to this result, Eq. (HHi, Fig. 1^ which represents Tj-fhAur also corresponds to the fre¬ 
quency dependence of the factor r on a logarithmic scale. 

Near ujqi the ratio r is very small (about a 0 . 25 %) and the secondary peak cannot be appreciated 
in the main part of Fig. |^a). In the inset we show the same data on a large vertical scale, in order to 
appreciate more clearly this secondary peak, whose temporal separation with the principal peak is 
about 355.36 ps. Near uo 2 the parameter T^/hAuj^ is lower than in the previous case and the ratio r 
is greater (roughly 40 %) and several peaks can be appreciated, as in Fig.j^b). It is clear that near 
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Figure 5: Transmitted pulses through our ID grating for an incident gaussian pulse of width 68.4 ps and a 
central frequency equal to cuoi = 1754 THz (a) and ojo 2 = 2225 THz (b). The inset in (a) corresponds to 
the same curve as in the main part with a different vertical scale. 


the transmission window edge the importanee of seeondary peaks is signifieant. Their separation, 
1297.80 ps, is larger than in the eenter of the window beeause this is inversely proportional to the 


frequeney distanee between two eonseeutive Lorentzian resonanees Acer (see Eq. (15)). 

In view of the previous results, for the high eontrasts eonsidered, it is eonvenient to use pulses 
with its eentral frequeney near the eenter of the transmission window (so as to reduee as mueh as 
possible the ratio between eonseeutive peaks), despite the faet that the group veloeity Vg is higher 
in this region. In Fig. we eompare the shape of a free pulse of width 6.84 ps (dashed line) 
with its eorresponding transmitted pulse (solid line). The eentral frequeney of the ineident pulse 
is ojqi = 1754 THz. This is an ultrashort temporal pulse (bandwidth of 413 GHz) desirable for 
teleeommunieation systems We ean observe a large delay of 12 pulses without signifieant 

pulse broadening and high transmission. We will return to this point in the next seetion. 
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Figure 6: Free propagation pulse of width 6.84 ps (dashed line) and its corresponding transmitted pulse 
(solid line) for a central frequency of woi = 1754 THz which traverses the ID grating described in Fig.[T] 

IV. BANDWIDTH DEPENDENCE OF THE TRANSMITTED PULSE PARAMETERS 

Let us now study numerically how significant pulse parameters, such as its transmission coef¬ 
ficient and its fractional delay, depend on the frequency bandwidth of the incident pulse 2Ace. 

In Fig. I^a) we represent the transmission coefficient |f 7 vP as a function of 2 Ace for an inci¬ 
dent pulse whose central frequency is cuoi = 1754 THz. The left dashed line corresponds to the 
frequency width of the Lorentzian resonance which is equal to 16.8 GHz (see Fig. |^a)) while the 
right dashed line (approximately 34 GHz) is the lower value of 2Aci; for which the transmission 
coefficient |fArp is independent of the frequency bandwidth. One observes that for short band- 
widths the transmission coefficient is practically unity because the frequency pulse is narrower 
enough to lay within a transmission resonance. As we increase the bandwidth 2Aa; the frequency 
pulse extends over several resonances what results in a lower transmission coefficient. When the 
number of resonances within the frequency pulse is large enough, the transmission coefficient 
Ifivp reaches a constant value independent of the frequency bandwidth, in this case 0.86. This is 
the regime where the traversal time is well approximated by the classical expression, Eq. (|^. For 
values of 2Aa; greater than 34 GHz the transmitted pulse consists of a serie of decreasing pulses 
as described in the previous section (see Fig.j^a)) with a short ratio r = 0.25%. 

Similar comments can be applied to Fig.|^b) where the central frequency of the incident pulse 
is now equal to uJo 2 = 2225 THz. The frequency width of the Lorentzian resonance (left dashed 











11 



Figure 7: Transmission coefficient |fArp as a function of the bandwidth 2Aw for an incident pulse whose 
central frequency is wqi = 1754 THz (a) and ojq 2 = 2225 THz (b). In both figures, the left dashed line 
corresponds to the frequency width of each Lorentzian resonance and the right dashed line is the lower value 
of 2 Aw for which the transmission coefficient is independent of the frequency bandwidth. 

line) is 0.7 GHz (see Fig. [^b)) while the bandwidth limit (right dashed line) is approximately 
6.6 GHz. One observes that the transmission coefficient |tArp reaches a lower value than in the 
previous case. The reason is that, as explained in Sec. II, the overlap between Lorentzians is 
negligible near the edge and the transmission in the midpoint between peaks reduces considerably. 


As regard the fractional delay, we show in Fig. [^this parameter versus 2 Aw for wqi (continuous 
line) and wo 2 (dashed line). We have performed these calculations via Eq. Q though similar results 
can be obtained with a classical treatment, Eq. Q. It can be appreciated that the fractional delay 
is, as expected, greater for values of Wq near the frequency window edge (where the group velocity 
is lower). However the distortion of the transmitted pulse, which is related to the factor r, is 
considerably high in this region (see Fig.|^b)). 

In view of the previous results one can conclude that the distortion of the transmitted pulse is 
unacceptable large when its central frequency lays near the edge of the transmission window. So, 
it is convenient to work near the center of this frequency window, even though it corresponds to 
larger group velocities. 
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Figure 8: Fractional delay versus the frequency bandwidth of the incident pulse 2Aw for two different 
values of the central frequency, wqi (continuous line) and a;o 2 (dashed line). In the latter case, the ratio r is 
quite large. 


V. DISCUSSION AND CONCLUSIONS 

We have analyzed the transmission properties of pulses through one-dimensional periodie strue- 
tures with alternating index of refraction ni and n 2 . We also explored the best conditions to achieve 
the maximum delay with the minimum possible pulse distortion. 

As we have already mentioned, in the present work we have not included either absorption or 
layer variability. We think that this is well justified for most situations of interest. In a typical de¬ 
vice employed in slow light experiments, the losses due to absorption are negligible up to distances 
of the order of 2 m ifTTI . As the structures considered here are 3cm long, we do not expect any 
significant effect due to absorption. As regard of layer variability, it is more difficult to estimate its 
effects, but the following argument convinced us that they again do not need to be considered in 
first instance. MacGum et al. [|22ll studied the localization length associated with layer variability 
and concluded that its value decreases (and the importance of its effects increases) as we approach 
the band edge. They performed a numerical simulation for a layer width variability of 20 % and 
obtained a localization length of roughly 20 times the layer period. Lhomme et al. ll^ estimated 
that the typical layer width variability in the type of devices considered by us to be about 0.2 %. 
As in one-dimensional systems the localization length depends on the inverse of the square of the 
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disorder, we estimate that the minimum loealization length will be approximately 2 ■ 10^ layer 
spaeings. The deviees eonsidered have 40000 periods and so we expeet a negligible deerease of 
the pulse amplitude due to loealization by disorder. Near the band edge loealization effeets are 
potentially more important, due to shorter loealization lengths, so one should be eareful with large 
layer variability in this region. 

For our ideal model, the transmission eoeffieient In can be well approximated by a sum of 
Lorentzian resonanees and the ratio between their width Fr and their separation AcUr is a signifieant 
parameter to eharaeterize the distortion of the transmitted pulse, whieh eonsists of an exponentially 
deereasing series of similar pulses given by Eq. ( [T5] ). This pulse distortion is greater near the edge 
of the transmission window, so it is eonvenient to work with pulse frequeneies near the eenter 
of this window. Our numerieal ealeulations show fraetional delays of 12 pulse widths, without 
signifieant distortion and high transmission, for 6.84 ps width gaussian pulses. 


Acknowledgments 

The authors would like to aeknowledge Vladimir Gasparian for many interesting diseussions. 
M.O. would like to aeknowledge finaneial support from the Spanish DGI, projeet FIS2006-11126. 

Appendix A: Transmitted pulse through a periodic arrangement of Lorentzian resonances 


Let us eonsider that the transmission amplitude of our strueture tp^{u) is a periodie funetion in 
u as given by Eq. ([T^. The transmitted wavefunetion 4/tr(f) will then be given by 


/ +00 

duj exp[-iujt] <l>inc(ce) ^ /(w - Auj^n, FJ, 

-OO ^ 


(Al) 


where G is a normalization eonstant. After introdueing the parameter v = uj — Auj^n and perform¬ 
ing some trivial ealeulations one ean write Eq. ( |A1[ ) as 

r+oo 


/ +00 

(iu/(u, rr)$inc('y + Auj-n) exp[—inf]. 

-OO 


(A2) 


Eor a large pulse over the frequeney domain, the wavefunetion <I)inc(cu) is praotieally eonstant 
(so we ean extraet this term out of the integral in the previous equation). We ean also express the 


first exponential in Eq. (A2) as an integral of the delta funetion and obtain 


= cy: 


r-l-oo 


du 6{u — Acoj-n) exp[—icef] 
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X 



dvf{v, Fr) exp[—ift] 


(A3) 


Taking into account the expression for the Fourier series of a Dirac comb 


^ (5(a; - Awrn) = ^ exp 

n ACCi- „ 


27riuj 


n 

AcCr 


(A4) 


and performing the inverse Fourier transform of the Lorentzian f{v,Tj. 


f + OO - 

/ dvf{v,Tj.)exp[—ivt] = y^ 

J — OO 


exp 


m 

2h 


(A5) 


we have for the transmitted wavefunetion 




r,|t| 

2h 


E 


(*+00 


du exp 


lU 


/27in 

\AUr 


<j). 

^ 1 


[OJ, 


(A6) 


Identifying the term within braekets in Eq. ( |A6| ) as the inverse Fourier transform of the incident 
wavefunetion \l/inc(t —STrn/AcCr) multiplied by 2n, we ean finally derive the expression for 
Eq.^. 

One observes that the transmitted wavefunetion \I/tr(f) is formed by a series of equally spaeed 
peaks, separated by the inverse of Aur, and modulated by an exponentially deeaying faetor 
exp[—Frit 1/2^]. The ratio between the intensity of two eonseeutive peaks r is then given by 
Eq. ([^. 
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